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Convergence conditions for quantum annealing arc derived for optimization problems rep- 
resented by the Ising model of a general form. Quantum fluctuations are introduced as a 
transverse field and/or transverse ferromagnetic interactions, and the time evolution follows 
the real-time Schrodinger equation. It is shown that the system stays arbitrarily close to 
the instantaneous ground state, finally reaching the target optimal state, if the strength of 
quantum fluctuations decreases sufficiently slowly, in particular inversely proportionally to 
the power of time in the asymptotic region. This is the same condition as the other im- 
plementations of quantum annealing, quantum Monte Carlo and Green's function Monte 
Carlo simulations, in spite of the essential difference in the type of dynamics. The method 
of analysis is an application of the adiabatic theorem in conjunction with an estimate of a 
lower bound of the energy gap based on the recently proposed idea of Somma et. al. for the 
analysis of classical simulated annealing using a classical-quantum correspondence. 
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Quantum annealing (QA) recently attracts much attention as a novel algorithm for op- 
timization problems. 1-4 A fictitious kinetic energy of quantum nature is introduced to the 
classical system which represents the cost function to be minimized. The resulting system 
searches the phase space by means of quantum transitions, which are gradually decreased as 
time proceeds. If the initial state is the ground state of the initial quantum Hamiltonian, the 
system is expected to keep track of the ground state of the instantaneous Hamiltonian under 
a slow decrease of quantum fluctuations. From this viewpoint, QA is also called quantum adi- 
abatic evolution. 5 Most of the numerical studies 1-3,6-15 showed that QA is more efficient in 
solving optimization problems than the well-known classical algorithm, simulated annealing 
(SA). 16 ' 17 

Convergence theorems for stochastic implementations of QA have been proved for the 
transverse- field Ising model. 18 A power- law decrease of the transverse field has been shown to 
be sufficient to guarantee convergence to the optimal state for generic optimization problems. 
This power-law annealing schedule is faster than that of the inverse-log law for SA given in 
the theorem of Geman and Geman. 17 ' 19 However, these theorems for QA were proved for 
stochastic processes to realize QA. It has been unknown so far what annealing schedule would 



I/O 



J. Phys. Soc. Jpn. Full Paper 

guarantee the convergence of QA following the real-time Schrodinger equation. We have solved 
this problem on the basis of the idea of Somma et. al. 20 These authors found that the inverse- 
log law condition for SA can be derived from the adiabatic theorem for a quantum system 
obtained from the original classical system through a classical-quantum mapping. Although 
they also discussed some aspects of QA, their interest was to use quantum mechanics to 
simulate finite-temperature classical statistical mechanics. We point out in the present article 
that the convergence condition of genuine quantum annealing, in which the system follows the 
real-time Schrodinger equation without temperature, can also be derived by a similar analysis. 

Let us suppose that the optimization problem one wants to solve can be represented as 
the ground-state search of an Ising model of general form 

N 

H pot = - X] Ji(T i ~ Ji 3 a i a j ~ J ijkViVj&k , (1) 

i=l ij ijk 

where of denotes the z component of the Pauli matrix at site i. Quantum annealing is realized 
typically by the addition of a time-dependent transverse field 

N 

i=i 

which may be regarded as the quantum kinetic energy to be compared with the potential 
energy (pQ). Initially the coefficient of the kinetic term T(t) is chosen to be very large, and the 
total Hamiltonian 

H(t) = H pot + H kin (t) (3) 

is dominated by the second kinetic term. The coefficient T(t) is then decreased gradually 
toward 0, leaving eventually only the potential term. Accordingly the state vector \ip(t)), 
which follows the real-time Schrodinger equation 

i±\m)=H(t)\m), (4) 

is expected to evolve from the trivial initial ground state of the transverse field ([2]) to finally the 
non-trivial ground state of eq. (pQ). The problem of central concern in the present paper is how 
slowly we should decrease T(t) to keep the state vector arbitrarily close to the instantaneous 
ground state of the total Hamiltonian ([3]), namely the adiabaticity condition, to achieve the 
goal of minimization of eq. (H|). 

The adiabatic theorem 21 provides the excitation probability at time t as 

dH(t) 



\(n(t)W))f 



n(t) 



dt 



o(t) 



(e n (t) - e (t)) 2 



(5) 



where \n(t)} is the nth instantaneous eigenstate of H(t) with the eigenvalue £ n (t)- We assume 
that |0(i)) is the ground state of H(t) and |<0(O)) = |0(0)). The probability © needs to be 
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arbitrarily small for the success of QA. We therefore evaluate an upper bound of 

dH(t) 



n(t) 



Ot 



0(f) 



(6) 

(e„(t) - eo(t)Y 

For this purpose we estimate the numerator and the denominator of eq. ©. As for the 
numerator, it is straightforward to see 



n(t) 



dH(t) 



Of 



0(f) 



< 



(7) 



since the time dependence of H{t) lies only in the kinetic term, which has A terms. Note that 
dr/di is negative. 

A lower bound on the denominator of eq. © can be evaluated using an inequality for 
a strictly positive operator. 22 First we recall that the Perron-Frobenius theorem states that 
a non- negative square matrix M has a real eigenvalue Ao satisfying |A| < Ao for any other 
eigenvalue A. If all the elements of M are strictly positive, M« > 0, its eigenvalues satisfy the 
stronger inequality, 22 



|A|< 



1 



K+ 1 



Ao, 



where k is defined by 



M ik 
max . 

i,j,k Mjk 



(8) 



(9) 



We apply the above inequality (jHJ) to the operator M = (£ max — H(t)) , where E ma _ x is 
the largest eigenvalue of i? po t- All the elements of the matrix M are strictly positive in the 
representation that diagonalizes of because E max — H[t) is non-negative and irreducible (that 
is, any state can be reached from any other state within A steps at most). In the asymptotic 
region t > 1 where T(t) <C 1, the minimum element of M, which is between two states having 
all spins in mutually opposite directions, is equal to N\T(t) N , where N\ comes from the ways 
of permutation to flip spins. Replacement of -ffkin by —N shows that the maximum matrix 
element of M has the upper bound (B max — -Emm + N) N , where -Emm is the lowest eigenvalue 
of i^pot- Thus we have 

(-Etaax — -Emin + N) 



K ~ NT{t) N 

Since the eigenvalue of H(i) is denoted by e n (t), eq. 



(£ n 



S n (t)) N < 



K+l 



is rewritten as 



eo(t)Y 



Substitution of eq. (fTU|) into the above inequality yields 

2( J Bmax - e (t))N\ 



e n (t)-e (t) > 



T(t) N = AT(t) 



N(E n 



Emin + A) 



N 



N 
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where we used k 3> 1 in the asymptotic time region where T(t) is very small. The coefficient 
A is estimated using the Stirling formula as 



2 V / 2^V(£max-£0(t)) ( N 



Ne N \E max -E min + Nj ' V ; 

which demonstrates that A is exponentially small for large N. 

Now, by the combination of the above estimates (JT]) and (|12p . we find that the sufficient 
condition for convergence is that the upper bound of eq. ([6]) 

(14) 

A 2 T(t) 2N dt { ' 

is arbitrarily small. By equating eq. (|14p to a small constant 5 and integrating the resulting 

differential equation, we find 

r(i) = (at)- 1 /(27V-i) ) (15) 

where a is exponentially small for large N and is proportional to 5. The transverse field should 
be decreased following this functional form or slower. Therefore the asymptotic power decay 
of the transverse field guarantees that the excitation probability is bounded by the arbitrarily 
small constant 5 2 at each time. 

The same discussions apply to quantum annealing using transverse ferromagnetic inter- 
actions in addition to a transverse field, 

H kin (t) = -f (t) (X>f + E of of ) • (16) 

\i=l ij J 

A recent study showed the effectiveness of this type of quantum kinetic energy. 23 A modifi- 
cation of the strictly positive operator to {E max — H(t)) N ' 2 in the above argument leads to a 
lower bound of the energy gap as a quantity proportional to T(t) N ' 2 . The resulting asymptotic 
annealing schedule is 

f^oct- 1 ^- 1 ), (17) 

which is faster than the case with the transverse field only, eq. (j 1 5 j) . This result implies 
that the additional non-zero off-diagonal elements of H(t) would widen the energy gap and 
accelerate the convergence of QA. 

In this paper, we have derived conditions for convergence of QA under the real-time 
Schrodinger dynamics using the adiabatic theorem. The asymptotic power-law annealing 
schedule (|15j) guarantees the adiabatic evolution during the annealing process at its final 
stage i > 1. This condition coincides with our previous results for stochastic implementations 
of QA. 18 It is remarkable that essentially different types of dynamics share the same condition 
for convergence. Note that the power decay derived above applies to the asymptotic region 
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t S> 1. At the initial stage, T(t) must be tuned following a different functional form to satisfy 
the adiabatic condition. 

For the adiabatic theorem to be applicable, the energy gap between the ground state 
and the first excitation state should be finite. The inequality (|12p implies that this condition 
is always satisfied in the transverse-field Ising model as long as the system size is finite. In 
the thermodynamic limit, of course, the gap may vanish at the critical point. We emphasize 
that the system size N is kept finite in the present paper because our purpose is to study 
optimization problems in which the number of elements is always finite. 
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Research Fellowships for Young Scientists. 
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